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Abstract. We present numerical relativistic magnetohydrodynamic and emis- 
sion simulations aimed to study the role played by the magnetic field in the dy- 
namics and emission of relativistic jets in Active Galactic Nuclei. We focus our 
analysis on the study of the emission from recoUimation shocks since they may 
provide an interpretation for the stationary components seen at parsec-scales in 
multiple sources. We show that the relative brightness of the knots associated 
with the recoUimation shocks decreases with increasing jet magnetization, sug- 
gesting that jets presenting stationary components may have a relatively weak 
magnetization, with magnetic fields of the order of equipartition or below. 



1. Introduction 

Thanks to the recent developments in the numerical modelling of magnetized 
relativistic jets it is now possible to study the influence of the magnetic field 
in the dynamics of relativistic jets. Computation of non-thermal (synchrotron) 
emission using relativistic magnetohydrodynamic (RMHD) simulations as in- 
puts allow us to obtain synthetic radio maps that can be directly compared 
with actual observations (see, e.g., Gomez 2002, and references therein). This 
comparison between simulations and observations is proving a powerful tool in 
the understanding of the physical processes taking place in the jets of Active 
Galactic Nuclei (AGN). 



2. eRMHD Simulations 



The simulations have been performed using a numerical code that solves the 
RMHD equations in conservative form and cylindrical coordinates with axial 
symmetry (see Leismann et al. 2005 for more details). Numerical fluxes between 
contiguous zones are computed using an approximate Riemann solver. Spatial 
second order accuracy is achieved by means of piecewise linear, monotonic re- 
construction of the fluid variables. Averaged values of the conserved variables 
are advanced in time by means of a Runge-Kutta algorithm of third order. The 
magnetic field configuration is kept divergence-free thanks to the implementa- 
tion of a constrained transport method. Using the RMHD results as inputs 
we compute synthetic radio continuum emission (total and linearly polarized) 
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maps that can be directly compared with observations. For this we integrate the 
transfer equation for the synchrotron radiation (for an optically thin frequency 
and in arbitrary units) taking into account the appropriate relativistic effects. 
See Gomez (2002), and references therein, for a complete description. 
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Figure 1. Panels show the logarithm of the pressure and the beam mass 
fraction (ratio of jet to external medium rest-mass densities) for Model A (a; 
top) and Model B (b; bottom). 



3. Setting up the numerical simulations 

The models arc characterized by the injection velocity into the external medium 
( i.e, Lorentz factor F), and the specific internal energy. Two more parameters, 
the beam to external medium ratios of pressure, K, and rest-mass density at 
the injection position, are used to fix the properties of the atmosphere through 
which the jet propagates. An ideal gas equation of state is used. We assume 
that the magnetic field has a helical structure (see Lind et al 1989; Komissarov 
1999), with axial, B^, and toroidal, B*^, components given by: 

( \f r <r- ) ^ f ^'^r^r/^rn) if r<r^ ] 

If r>r- ' Biirm/r) if r„, < r < r, (1) 

^ ^ { if r > J 

where rj is the jet radius, and is the radius at which the toroidal magnetic 
field reaches its maximum value B^. The magnetic field is characterized by the 
averaged pitch angle (j), which is defined as tancf) = B'f'/B^, and the magnetiza- 
tion parameter /3, which is the ratio of magnetic field to gas pressure. 
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4. The Influence of the Jet Magnetization 

In order to study the influence of the magnetization in the dynamics and emission 
we have performed simulations of overpressured {K=2) relativistic (F = 4) jets 
carrying a helical magnetic field with pitch angle <p = 65° and two different 
magnetizations at the jet inlet: Model A with P=0.1, and Model B with /3=3. 




Figure 2. Toroidal magnetic field profiles at different radius along and 
across the jet for Model A (a; top) and Model B (b; bottom). 



The dynamics of the jet is dominated by the tendency of the beam fluid to 
achieve pressure equilibrium with the external medium, causing overexpansions 
and overcontractions past equilibrium of the fluid leading to the formation of a 
pattern of recollimation shocks (see Fig. 1). We can study the role played by 
the magnetic field in the dynamics by analyzing the influence of the toroidal and 
axial components separately. The toroidal component (see Fig. 2) is relatively 
weak close to the jet axis but it can be significant close to its maximum strength. 
This component generates a force against the expansions and contractions of the 
fluid, which will be maxiiniiiii at radius r^. Studying the position of r„j in the 
transversal profiles at different radius along the jet (see Fig. 2) we see that it 
changes during the expansions and contractions of the jet, but remains close to 
the original position for Model A while for Model B moves to an inner radius 
with distance along the jet. The contribution of the axial component, which has 
a constant profile inside the jet, is to generate a magnetic tension that opposes to 
the transversal movement of the jet. Therefore both components contribute to 
the confinement, collimation, and stability of the jet. Hence a more magnetized 
jet will not expand initially as much as a jet with a lower magnetization, and 
the subsequent contraction will take place before in the grid length, resulting in 
weaker shocks. Certainly, by comparing our models A and B (see Figs. 1 and 
2) we observe that i) the recollimation shocks are weaker, and there is larger 
number of them, as f3 increases; ii) Model B presents a stronger collimation, 
induced by the enhanced magnetic field. 
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Figure 3 shows the computed synthetic emission maps for a viewing angle 
of 14° for models A and B. We observe bright knots of emission associated with 
the recollimation shocks, which may provide an interpretation for the stationary 
components commonly seen in multiple sources (e.g., Gomez et al. 1995; Agudo 
et al. 2001; Jorstad et al. 2005). Figure 3 shows that the relative brightness of 
the knots associated with the internal shocks decreases with increasing magneti- 
zation. Indeed, we observe that Model A presents strong stationary components 
with a relative brightness that doubles the emission of the underlying flow. On 
the other hand. Model B, which has very weak recollimation shocks (see Fig. 1), 
shows a smooth jet emission with small variations along the jet. We can there- 
fore suggest that jets presenting stationary components may have a relatively 
weak magnetization, with (3 close or below equipartition. 



a) b) 




Figure 3. Panels show total emission maps in arbitrary units (top) and total 
emission along the jet axis direction (bottom) for a) Model A and b) Model 
B for a viewing angle of 14°. 
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